An analysis of non-classical austenite-martensite interfaces in CuAINi 
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Abstract 

Ball & Carstensen [31 [3], theoretically investigated the 
possibility of the occurrence of non-classical austenite- 
martensite interfaces and studied the cubic-to-tetragonal 
case extensively. Here, we aim to present an analy- 
sis of such interfaces recently observed by Seiner et al. 
[12"] in CuAINi single crystals, undergoing a cubic-to- 
orthorhombic transition. We show that they can be de- 
scribed by the non-linear elasticity model for martensitic 
transformations and we make some predictions regarding 
the volume fractions of the martensitic variants involved, 
as well as the habit plane normals. 

Introduction 

A classical austenite-martensite interface is one in which 
the undistorted austenite meets a simple laminate of 
martensite and these have been broadly studied. In recent 
years, a theory of martensitic transformations has been 
developed which allows austenite-martensite interfaces to 
occur in which the microstructure of the martensite is 
more complicated; these are referred to as non-classical 
interfaces and have hitherto not been systematically ob- 
served. In this non-linear elasticity model, in which inter- 
facial energy is neglected, microstructures are identified 
as limits of minimizing sequences for the total free energy 

Ie{y)= [ <p(Vy(x),6)dx. (1) 
Jn 

Here, f2 represents the reference configuration of undis- 
torted austenite at the transformation temperature 9 C and 
y{x) denotes the deformed position of particle x 6 f2. 
The free-energy function <p(F,0) depends on the defor- 
mation gradient F and the temperature 9. By frame 
indifference, <p{RF, 9) — tp(F, 9) for all F, 8 and for all 
rotations R; that is for all 3 x 3 matrices in SO(3) = 
{R : R T R = l,dcti? — l}. By adding an appropriate 
function of 9, we may assume that minp ip(F, 9) — 0. At 
9 C , the energy wells of the free-energy function are given 



by SO(3) for the undistorted austenite and SO(3)Ui, 
i = 1,...,N, for the N distinct variants of martensite, 
where each Ui is a positive definite, symmetric matrix. 
Hence, <p(F, 9 C ) > with equality if and only if F E 
SO{3) U U£Li SO(3)Ui. For 9 < 9 C , the martensite wells, 
Ui=i SO(3)Ui, minimize <p, whereas for 9 > 9 C the min- 
imum is given by the austenite well, SO(3)a(0)l, where 
a{9) is the thermal expansion coefficient of the austenite 
and a(9 c ) = 1. 

A non-classical planar austenite-martensite interface 
{x ■ to — k} corresponds to a choice of habit-plane nor- 
mal to such that there exists a sequence of deformations 
yW for which I$ c {% ' ) — > as j — > oo. We require 
that, for x • to < k, the values of the deformation gradient 
Vy^' tend to SO(3), i.e. yU> corresponds to the austen- 
ite phase; without loss of generality, this is equivalent to 
Vi/-'' (x) — 5- 1 except possibly for a set of zero volume. 
For x-m > k, we require that, as j — > oo, Vj/- 3 ' tends in 
a suitable way to the set K = U i=1 SO(3)Ui of martensite 
energy wells. In fact, we require that the Young measure 
{v x )xen of is supported in the set K (for details see 

e.g. P). 

From now on, we shall make the assumption that the 
martensitic microstructure is homogeneous; that is, for 
x ■ m > k, the macroscopic deformation gradient F = 
Vy(x) is independent of x. The set of all such matrices 
F is called the quasiconvexification of K and is denoted 
by K qc . It can be shown that K qc is also the set of F 
such that there exists a sequence of deformations z^' with 
z^'(x) — Fx on the boundary of il, dfl : and V z^\x) — > 
K in the above sense. 

If we know K qc for a given set of martensite wells, to en- 
sure geometric compatibility, we need to examine whether 
it is possible to establish a rank-one connection between 
SO(3) and K qc ; that is, we need to find vectors b and m 
such that 

1 + 6® to G K qc , (2) 

where, by frame indifference, we have chosen the iden- 
tity matrix 1 to represent the austenite energy well. 



Unfortunately, we only have a characterization of K qc 
for two martensitic wells (N — 2), i.e. when K — 
S0(3)U! U SO(3)U 2 . In this case, any F G K qc can be 
obtained as the macroscopic deformation gradient of a 
double laminate (see [S]). Even though K qc is unknown 
in the case of three tetragonal wells, Ball & Carstensen 
[3] were able, using the two- well calculation, to character- 
ize exactly the values of the deformation parameters for 
a cubic-to-tetragonal transformation which permit non- 
classical interfaces. They also presented results for a 
cubic-to-orthorhombic transformation, though not of the 
type occurring in CuAlNi. 

Given a matrix F, the question as to whether we can 
solve the equation 

1 + 6® to € SO(3)F (3) 

for vectors b and to is answered by the following lemma 
in [H[8]. 

Lemma 1. Let F be a non-singular matrix that is not 
a rotation. Then the wells SO(3) and SO(3)F are rank- 
one connected if and only if the middle eigenvalue of F T F 
is 1. Then 1 + b ® m £ SO(3)F for some b if and only if 
m is a non- vanishing multiple of one of the two vectors, 
\/l — Aiei ± — le3, where < Ai < 1 < A3 are the 
three eigenvalues of F T F with corresponding eigenvectors 
ei, e 2 , e 3 . 

Having outlined a brief description of the model, we 
proceed to the experimental observations on a CuAlNi 
single crystal. 

Experimental Observations 

The first micrographs of interfaces between austenite and 
crossing twins were obtained by Seiner et al. [13] , who 
documented that such interfaces can form during the 
shape recovery process of CuAlNi single crystals. How- 
ever, these interfaces were only weakly non-classical, i.e. 
they were classical interfaces weakly disturbed by a negli- 
gible volume fraction of compound twins intersecting the 
first-order laminate of the Type-II twins. These obser- 
vations motivated the authors towards more intensive re- 
search in this field - the experiment was improved in order 
to increase the volume fraction of the compound twins in 
the martensitic microstructure. The resulting experimen- 
tal procedure is briefly outlined below, and will appear in 
more details in [T2"] . 

The specimen examined in this case was a 
3. 9x3. 8x4. 2mm rectangular parallelepiped of the 
austenitic phase of CuAlNi, cut from a single crystal 
of this alloy such that the normals to the specimen 
faces had approximately the principal crystallographic 



directions (100). The original single crystal was grown 
by a Bridgman method at the Institute of Physics ASCR 
in Prague. The specimen was subjected to the following 
sequence of mechanical and thermal loadings (see Fig. IT]) : 



COMPRESSION 




Figure 1: Outline of the experimental procedure. 

(a) At room temperature, the specimen was transformed 
into a single variant of 2H martensite by applying 
uniaxial compression (in a bench vice). Due to the 
effect called mechanical stabilization of martensite, 
the specimen did not return to austenite after un- 
loading, but remained as a single variant of marten- 
site. 

(b) The specimen was rotated by 90° and uniaxial com- 
pression was applied again. In this case, the load- 
ing induced the reorientation of martensite into an- 
other variant via compound twinning (for an anal- 
ysis of the relation between the direction of applied 
compression and the activated twinning systems in 
CuAlNi, see [§])■ The reorientation was not fully 
completed. Instead, the loading was interrupted at 
the moment when the specimen contained compa- 
rable volume fractions of both variants. By such a 
procedure, we obtained a finely compound twinned 
specimen. 

(c) The finely compound twinned specimen was heated 
from one side using a gas lighter, which induced the 



shape recovery process, i.e. the thermally driven re- 
turn of the specimen into austenite. As the com- 
pound twins cannot form any compatible interface 
with austenite, the transition was achieved by for- 
mation of an interfacial micro structure, which en- 
sured a compatible connection between the mechan- 
ically stabilized martensite (the compound twins) 
and austenite. This interfacial microstructure was 
formed by Type-II twins crossing the original com- 
pound microstructure and getting arranged into a so 
called X-interface (for more details of formation of X- 
interfaces in CuAINi see [13 , for the theoretical anal- 
ysis of this microstructure, see [10]). By removing 
the heating in the middle of the course of the transi- 
tion, the interfacial microstructure was stopped, and 
the non-classical interfaces between austenite and the 
two mutually crossing systems of twins (compound 
and Type-II) were observable by optical microscopy. 
An example is given in Fig. [2j 



orientation of the habit plane varies as well, as will be 
shown in the theoretical analysis given in the following 
section. However, a more complete analysis of the curved 
interface remains to be done. 




Figure 2: Optical micrograph of a non-classical inter- 
face between austenite and a martensitic microstructure. 
The arrows indicate the orientations of twinning planes 
of Type-II and compound twinning systems. 

This procedure was repeated several times, which enabled 
the capturing of several optical micrographs of the non- 
classical interfaces. An interesting observation was that, 
as shown in Fig.[3j the interface between the austenite and 
the system of crossing twins was never exactly planar, but 
rather slightly curved. This results from the fact that the 
pattern of compound twins induced in the specimen by 
the uniaxial compression in stage (b) of the experimental 
procedure is never exactly homogeneous. With varying 
volume fraction of the compound twins, the admissible 




Figure 3: Curved interface between crossing twins and 
austenite resulting from the inhomogeneity of compound 
twinning. (Optical microscopy.) 



Analysis of Non- Classical Interfaces 

In this section, we present an analysis for the above non- 
classical austenite-martensite interface via the non-linear 
elasticity model for martensitic transformations. The 
martensitic region on one side of the observed interface 
consists of twin crossings involving four martensitic vari- 
ants. Since the quasiconvexification of more than two 
wells is unknown, we are not able to analyze all the pos- 
sibilities for non-classical interfaces in CuAINi. However, 
we can do this in the twin crossing case. Firstly, we shall 
concentrate on the martensite phase and try to construct 
the microstructure consisting of compound and Type-II 
twin crossings, as shown in Fig. [4] 

Let Ua and Ub be two martensitic variants able to 
form a Type-II twin and Ua>, Ub> the respective com- 
pound counterparts. Clearly Ua' and Ub' can also form 
a Type-II twin and we proceed by writing the compatibil- 
ity equations for the parallelogram microstructure. These 
are: 

RabUb ~Ua = oab O n A B (4) 



1 nAB 

Ra'B'Ub' - U A > = b A 'B> ® n A 'B> 
Raa'Ua' -U a = b A A> <& n A A> 
Rbb'Ub' -U B = b B B> ® ubb' 



(5) 
(6) 
(7) 




Figure 4: Parallelogram micro-structure 



where bjj and n/j are, respectively, the shearing vector 
and the normal to the twinning plane for the system of 
variants Uj and Uj. Also, Rjj denotes the mutual rota- 
tion between variants Ui and Uj. 

If the above compatibility equations hold, a necessary 
and sufficient condition for the entire parallelogram mi- 
crostructure to be compatible is that 



RabRbb 1 — Raa'Ra' 



(8) 



Equation ([8| is necessary, since both sides describe the 
mutual rotation between U A and Ub> , and sufficiency fol- 
lows by showing that Q-([7| imply that the normals uab, 
nA'B'i tiaa' and tibb' are coplanar. 
Let 

M AB = (1- X)U a + XRabUb (9) 



Ma'B' = (1 - X)U A > + \Ra-b-Ub- 



(10) 



which represent the macroscopic deformation gradients 
corresponding to the Type-II structures. From [7J, there 
are solutions of the twinning relations Q and Q with 



riAA> = n B B> 



(11) 



and we make this choice in accordance with the exper- 
imental observations. Then the geometry of the paral- 
lelogram structure requires A to be the same in ^ and 

At this stage we note that det M A > b 1 = det Mab ■ To 
see this, consider equation Q. We get that 

RabUb = Ua{1 + U^bAB <S) u A b)- 



Therefore, by taking determinants, we obtain 

det U B = det U A {1 + U^^ab ■ tiab)- 
Clearly det Ub = det Ua and thus 

U^bAB ■ n AB - 0. (12) 
Hence, for Mab, we have 

det Mab = dctlf A - 

Similarly, det Ma> b' = det Ua> = det Ua and we deduce 
that det M A ' B' = det Mab ■ This result has a simple 
physical interpretation. Since the determinants have the 
meaning of volume change and the microstructures Mab 
and Ma> b 1 are just mixtures of variously rotated single 
variants, the volume must remain the same. 

Now, assume that the volume fraction of variant U a' in 
Ua, as well as that of Ub' in Ub, is A. Then, the macro- 
scopic deformation gradient of the entire microstructure 
is 

M = (1 - A) Mab + KR A A'M a >b-- (13) 

Whether the interface between the parallelogram mi- 
crostructure and the austenite can be formed is a ques- 
tion of finding A and A such that the middle eigenvalue 
of M T M is equal to 1, which is equivalent to finding a 
rank-one connection between the austenite and SO(3)M. 
In particular, by Lemma 1, 



det(M T M 



0. 



(14) 



Forming the Interface 



By manipulating the compatibility equations and making 



use of equations (JSj) and (11) we can see that 

(bAA>-RABbBB>)®nAA'-bAB®nAB + RAA'bA'B>®nA>B' 

= 0. 

All normals cannot be parallel to each other and hence, 
using a result in [10], we deduce that 

bAA' - RABbBB' || b A B \\ RAA'bA'B'- (15) 

Hence, there exists some r\ such that 

J]b A B = b A A' - RABbBB'- (16) 

By identifying the rotation Rab and using the formulae 
for the twinning solutions in Result 5.2 of [5], as well as 
the relations between them [7J , we obtain that 



V = 



\b 



-,bAA' ■ b A i 



AB 



Our goal is to deduce an expression for det (M T M — 1) 
which will enable us to find solutions of ( |l~4"| for A, A £ 
(0, 1). Towards this end, we shall seek a rank-one con- 
nection between M A b and Ma'b 1 - Indeed, using ([8| and 
equations Q-([7]), 

RaA'Ma'B'-M A b = {l-\)b A A'®nAA'+^RABbBB'®n B B' 

(17) 



-AU A 1 n A A' ® U A 1 b A A> 
+A\UA 1 n A A> ® E/^&ab^a^ab • + ??). 
Here we have already used the fact that 

U^UAA' ■ b A B = 0. 



(22) 



Having that %aa' — kbb' and using (16 1 we get 

Raa'Ma'B' - M A b = b 
where 6* = bAA' — Xr/bAB 



TlAA' 



Combining equation (17) with (13) we get that 



M = M AB + A6* ® n A A> 
and using the expression for Mab we obtain 

M = U A + Xb A B ® n A B + A&* ® n^A'- 

Then, 

dct(M T M - 1) = det M T det(M - AT T ). 
We shall first calculate det M T = det M. Hence we have, 
det M = det Mab det(l + AM AB b* (8 tua')- 



(18) 
(19) 



However, from (17) 



Raa'Ma-b' = M AB (1 + M AB b* ® tua') 
and taking determinants, we see that 

det M A 'B' = det Mab(1 + M AB b, ■ n AA >)- 
Since det M A > B' = det Mab we get that 
Mab & * ■ n A A> = 0. 



(20) 



From (19), we have 



det(M T M - 1) = det U A det(M - AT 



and it remains to calculate M — M 



-T 



Making use of the expression for M and ( 20 1 



M~ T = M- A l - AM^n AA ' <8 M AB b, 



and from ^ and Q 

m ab = Ua 1 ~ W^bAB ® C^ab- 

Combining the last two equations and after some calcu- 
lations wc obtain that 



M~ T = U^ 1 



XU A l n AB ®U A 1 b 



A U AB 



(21) 



The above result requires some effort and is provided by 
investigating the axes of the rotations in the austenitic 
point group along with Result 5.2 in |B]. 

Using the expressions for M and M~ T we get 



det(M T Af - 1) = dct(A + XAi + AA 2 + XAA 3 
where 



(23) 



A 
At 

A 3 



Ui 



1 



(24) 
(25) 
(26) 



UAbAB ® riAB + n A B ® C^4 b A B 
U A bAA> ® TlAA' + n A A> <& t^4 1 !*AA' 

-(Cl^xa • &AA' + ifKtA' ® t^&AB (27) 

-rjU A bAB ® n AA' • 



It is trivial to observe that, for fixed A, the expression 
multiplying A is of rank 2 and similarly, for fixed A, the 
expression multiplying A is of the same rank. Hence, the 
determinant of M — M~ T will only contain terms with A 
and A in powers not greater than two. Letting g(X, A) := 
det(M T M - 1) we deduce that 

g(X, A) = a + (3X + 7 A + aX 2 + bA 2 + cXA + dXA 2 (28) 

+eA 2 A + /A 2 A 2 . 

Using the relations between the twinning solutions for dif- 
ferent martensitic variants [7J, the above expression sim- 
plifies further by noting that for all A £ (0, 1) 



«7(O,A)=0(1,A) 
and, respectively, for all A £ (0, 1) 

3 (A,0)= 5 (A,1). 



(29) 



(30) 



Now, from equations (|28j)-([30j) we obtain the following 
form for g; 

g(X,A) = a + ai (A 2 -A) + a 2 (A 2 -A)+a3(A 2 -A)(A 2 -A). 

(31) 

From (23) and (31) it becomes easy to specify the coeffi- 



cients ao,ai,<22,a3 in the expression for g(X, A). Firstly, 
a =5(0,0)=det(C7i-l). (32) 



Noting that a\ = — gf(0, 0) and a 2 - 
duce, after a simple calculation, that 



OA 



(0,0), we de- 



ai = -2b AB ■ U A coi{U\ - l)n AB 



(33) 



equal to 1. This is answered in a similar fashion to @] 
and hence, we shall require that 



trUl 



-det Ul- 



-2 + (A 2 -A) \b AB \ 2 + (A 2 -A) \b AA ,\ 2 (41) 



and 



a 2 



> A - 1)«AA 

expression gets more complicated reducing to 



(34) 



2b AA , ■ U A coi(U 2 A 
As for the last coefficient, a 3 = 4(||(0, |) + ai) and the 



a 3 = 4cof(A + X -A 2 ) ■ {A x + ~A 3 



4ai. 



(35) 



For A = 0, we only have the Type-II structure and we get 
g(X,A) = a a + ai {\ 2 - A). 



This becomes zero for A 2 — A = - M ; call this A = A* e 
(0, |] , which agrees with the value obtained, for example, 
by Ball & James in pQ. 

Setting g(A, A) = and solving for A we deduce that, 
as long as 



then 



Ol 

A 2 — A = - 



[0, 



q +a 2 (A 2 - A) 
oi +« 3 (A 2 - A)' 



(36) 



(37) 



For A = and A = 1 the above equation admits two 
solutions, namely A* and 1 — A*, and we see that branches 
starting from (A*,0) and (A*,l) are created consisting 
of values of (A, A) that make the non-classical interface 
possible. Provided that 

ao«3 7^ a ia>2 (38) 
a necessary and sufficient condition for these to meet, i.e. 



for equation (37) to have solutions for all A g [0,1], is 
that 

4an — &9 1 , . 

< -° < -. (39) 

4ai — a 3 4 



If, in addition, 



4a n 



4a i 



«2 , 1 
03 4 



(40) 



+ (A 2 -A)(A 2 -A)r ? 2 |6 ylB | 2 >0 

for all pairs (A, A) that make the interface possible. This 
will also imply that the other eigenvalues of M T M, Ai 
and A3, are bounded away from 1, i.e. < Ai < 1 < A3. 

However, the formulae for the cubic-to-orthorhombic 
transformation get too involved and for this reason we 
will proceed numerically. 

Numerical Results 

In the remainder of this section we present a numeri- 
cal calculation where, in accordance with [TT] , the lattice 
parameters for CuAINi were chosen to be a = 1.06372, 
(3 = 0.91542 and 7 = 1.02368. The martensitic variants 
used are the ones obtained from the experimental obser- 
vations; that is, A = 3, B = 6 with A' = 4, B' — 5 their 
compound counterparts. Following the above analysis, 
we calculated the zeros of g, i.e. the values of (A, A) that 
allow the interface to occur. Relations (36) and (38)-(41) 



were satisfied and we plotted A against A as shown in 
Fig-i 
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Figure 5: Values of A that make the interface possible for 
< A < 1 



then the solutions will be distinct. 

Otherwise, if a^a^ = a.\ot2, g simplifies even further 
and we obtain that, for all A g [0,1], A = A* and A — 
1 — A* will suffice. Due to the symmetry of g, branches 
are also created at the points (1 — A*,0) and (1 — A*,l) and 
the condition for these to meet is the same. A remaining 
question is whether it is the middle eigenvalue that is 



It is easily seen that the value of A does not change sig- 
nificantly from A* (corresponding to A = 0), which would 
give the classical interface between the Type-II twinning 
system and the austenite. 

Moreover, using the algebraic procedure given in [4] , we 
calculated the different normals m(A, A) for (A, A) on the 
curves of Fig. [5] and a plot of these is given in Fig. [6j 



where the normals are depicted as points on the unit 
sphere. These normals lie on four segments of curves 



-1.0 




between this volume fraction A and a volume fraction of 
the Type-II laminate A intersecting compatibly the com- 
pound twins was derived analytically for a general cubic- 
to-orthorhombic transition, the analysis brought by this 
paper can be easily applied to predict the existence of 
similar non-classical interfaces in any other shape mem- 
ory alloy with the same class of transition (e.g. CuAIMn) 
as well as for materials undergoing the cubic-to-tetragonal 
transitions (since the tetragonal symmetry is a member 
of the orthorhombic symmetry class.) 

The numerical simulations carried out in the last sub- 
section of this paper reveal that there is a dependence be- 
tween the compound volume fraction and the habit plane 
orientation for the CuAINi alloy. This finding is consis- 
tent with the optical observations of slightly curved non- 
classical interfaces between austenite and the crossing- 
twins microstructure with heterogeneous compound vol- 
ume fraction (Fig. [3|. 
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Figure 6: Plot of the normals on the unit sphere for (A, A) 
as in Fig. [5] 
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Conclusion 

This paper brings a theoretical analysis of compatible 
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tion of the compound laminate. Since the relation (37) 
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